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Nonequilibrium condensate systems such as exciton-polariton condensates are capable of support- 
ing a spontaneous vortex nucleation. The spatial inhomogeneity of pumping field or/and disordered 
potential creates velocity flow fields that may become unstable to vortex formation. This letter 
considers ways in which turbulent states of interacting vortices can be created. It is shown that 
by combining just two pumping intensities it is possible to create a superfluid turbulence state of 
well-separated vortices, a strong turbulence state of de-structured vortices, or a weak turbulence 
state in which all coherence of the field is lost and motion is driven by weakly interacting dispersive 
waves. The decay of turbulence can be obtained by replacing an inhomogeneous pumping by a 
uniform one. We show that both in quasi-equilibrium and during the turbulence decay there exists 
an inertial range dominated by four-wave interactions of acoustic waves. 
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Introduction. The phenomenon of turbulence - chaotic 
motion of vortices of many different length scales - is 
ubiquitous in nature, and quantitative understanding of 
it is a notoriously difficult problem of classical physics. 
Turbulence occurs in many usual fluid flows as well as 
in exotic systems such as plasmas and superfluids. Vor- 
ticity in superfluids is quantized in units of h/m, where 
m is the mass of the boson in contrast with continu- 
ously distributed vorticity of a classical Navier-Stokes 
fluid. In superfluids quantized vorticity is considered to 
be an evidence for a macroscopically occupied quantum 
state that can be described by a classical complex-valued 
wave function ^/>(x, i). Quantization of velocity circula- 
tion in superfluids leads to significant differences between 
superfluid turbulence (ST) and classical turbulence. On 
the other hand, at large Reynolds numbers the motion 
of well-separated vortices in an incompressible classical 
flow may have similar features to ST. In this case the 
vortex dynamics in superfluids is almost classical in ac- 
cordance with the Biot-Savart law (BSL). The decay of 
the turbulence (loss of the vortex line density) occurs 
due to dissipative effects induced by interactions with a 
normal fluid component (with a thermal cloud). 

Recently by introducing an external oscillatory per- 
turbation in a trapped atomic BEC it became possible 
to obtain a disordered system of many topological de- 
fects pQ. The dynamics of this matter field differs from 
both dynamics of vortices in classical turbulence and in 
superfluid helium turbulence. Firstly, the characteris- 
tic distance between vortices is comparable to their core 
sizes, so the chaotic behavior is seen on the level of a sin- 
gle vortex, secondly, these vortices are not structured, so 
they do not obey BSL, finally, the system is in a strongly 
non-equilibrium state. These creates a novel nontrivial 
regime of a classical complex matter field — "strong tur- 
bulence" state - whose evolution is quite different from 
that of ordered condensate. In analogy with other nonlin- 



ear systems such as plasmas, fluids and nonlinear optics, 
apart from the regime of strong turbulence there exists 
the regime of weak turbulence where all phases of the 
complex amplitudes of the matter field are random. Re- 
cently |2 these three regimes (superfluid, strong turbu- 
lence and weak turbulence) have been observed at differ- 
ent temperatures in 2D cold atomic gases, showing a uni- 
versal scaling. The weak turbulence plays crucial role in 
kinetics of Bose-Einstein condensation [3] . It was shown 
that a strongly non-equilibrium Bose gas evolves from the 
regime of weak turbulence to superfluid turbulence, via 
states of strong turbulence in the long- wavelength region 
of energy space. An important question remains whether 
it is possible to force a condensate system to pass through 
these stages in a reverse order. It has been suggested [I] 
that if a sufficiently strong external perturbation is ap- 
plied to the trap, it is in principle possible to obtain the 
weak turbulence state. When this is done it will lead to 
a discovery of nontrivial transitional regimes of classical 
matter fields in atomic systems [5]. 



In the last few years the Bosc-Einstcin condensation 
has been achieved in solid state systems [6], such as mi- 
crocavities, ferromagnetic insulators and within super- 
fluid phases of 3 He. Microcavity exciton-polaritons are 
quasi-particles that consist of superpositions of photons 
in semiconductor microcavities and excitons in quan- 
tum wells. The Bragg reflectors confining photon com- 
ponent are imperfect, so exciton-polariton have finite 
life time and and have to be continuously re-populated. 
Such combination of pumping and decay leads to quasi- 
particle flow even at steady states of the system. At 
sufficiently low densities these quasi-particles can form 
a Bose-Einstein condensate, so the many particles quan- 
tum system can be described by a classical equation in a 
form of the complex Ginzburg-Landau equation (cGLE) 
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FIG. 1: (color online) The time snapshots of the density \i/j\ 2 
of the fields ip obtained by numerical integration of Eq. |T]) 
for 77 = 0, a = 0.3 and (i) a(x) = 2 for )x - (±5, 0)| < 2 and 
1/2 otherwise, t=100 (left panel) and (ii) V ex t = x 2 + y 2 and 
a(x) = 5 for x 2 + 2y 2 < 64 and a — —0.5 otherwise (right 
panel). Red ellipse indicates the pumping spot. Luminosity 
of the density plots is proportional to density. Vortices are 
seen as black dots. 
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2id t ^ = [-V 2 + + 1 (a - ir)d t Tp - o\i>\ 2 )\ (1) 

where a is an effective gain that represents intensity of 
the pumping field, a represents nonlinear losses. The 
unit of length is a healing length £ = h/y/2-mUpoo that 
defines the size of the vortex core and the unit of time 
is U/21/poo, where U is the strength of a <5-function in- 
teraction potential. We shall assume that a = ao is 
constant away from some localised nonuniformities and 
so the number density there is p^. It is possible to in- 
clude a disorder potential of the microcavity by adding 
14 x t(x)?/> to the right-hand side of Eq. Q. 

This equation is a mean-field description of the conden- 
sate; it can also be derived from the saddle point in a path 
integral formalism 10J. In the absence of pumping and 
dissipation Eq.Q reduces to the Gross-Pitaevskii equa- 
tion describing an equilibrium Bose-Einstein condensate. 
The energy relaxation has been noted to be of impor- 
tance in experiments on extended ID waveguides |HU12j . 
These effects can be included in the cGLE by means on 
a parameter 77 JT3] . This is the same term that has been 
incorporated into the Gross-Pitaevskii equation to rep- 
resent a dissipation of the condensate component due 
to interactions with a thermal cloud [T3]. The turbu- 
lence and mechanisms of the vortex generation in equi- 
librium condensates are well known. These include (i) 
interactions of finite amplitude sound waves (e.g. energy 
exchange between rarefaction pulses may lead to vortex 
formation) [T5]; (ii) existence of critical velocities of the 
flow (e.g. moving objects generate vortices if the Lan- 
dau critical velocity is reached on their surfaces [16]): 
(iii) modulational instabilities of density variations (e.g. 
transverse instability of a dark soliton in 2D generates 
vortices) [T7]. Some of these mechanism may produce 
vortices in the cGLE as well. For instance, the flow of 



exciton-polaritons about a spatially extended defect may 
produce vortex pairs of opposite circulation depending 
on the flow velocity [18] . In addition, Eq. ([I]) with gain 
and dissipation can form vortices by other physical mech- 
anisms. For instance, an inhomogeneity of the pumping 
or/and disorder potentials form steady currents which 
may produce vortices through pattern forming symme- 
try breaking mechanism [8j. 

Although the formation of vortices has been observed 
in experiments [7 they seem to appear due to the in- 
trinsic disorder potential in CdTe. The vortices become 
pinned at the local minimum of such potential and re- 
main stationary. So the conditions in which a turbulent 
state of matter can be obtained in exciton-polariton con- 
densates remained unclear. The purpose of this letter 
is to suggest how the turbulent state can be created in 
such a system, to study the properties and structure of 
the turbulence, and to propose how different regimes can 
be detected experimentally. It will be shown that turbu- 
lence can be created by deliberately designed pumping 
fields, and depending on characteristics of such fields the 
system can reach various regimes of turbulence from su- 
perfluid turbulence to strong and finally weak-turbulent 
state. 

Vortex formation. To illustrate the basic mechanism 
that drives the formation of vortices we first consider 
a pumping field in a form of a step function in ID, 
so that a = a\ + ao,a = o~\,r] = rji for x < and 
a = ao,a = cro,rj = t]q for x > 0. The steady state 
mass continuity and Bernoulli equations resulting from 
the Madelung transformation ip — ^fpexpiS applied to 
Eq. |l]) are p = u 2 + p — d 2 ^fp/l^fpdx 2 and d(pu)/dx — 
(a—r]p—ap)p where p is the number density, u = S'(x) is 
the velocity and the chemical potential p is introduced by 
2idtip = pip- Away from large density fluctuations we can 
drop the quantum pressure term d 2 - s /p/2 y /pdx 2 . We ex- 
pect that u — > as x — > —00, so p — > (on +a )/(cri +r/i). 
As x — > 00, therefore, there will be a steady current 
u= [a x (rj + cr ) + a (r/o - iji + a - G\ ) / cr (m + °i )] ~ 1/2 
generated by the step. The presence of boundaries or 
other sources of outflow generate interference fringes 
seen, for instance, in recent experiments in ID |llj . In 
2D the fringes that meet at nonzero angles evolve into a 
pair of vortices of opposite circulation as seen on the left 
panel of Fig. [T] The mechanism leading to vortex forma- 
tion in this case is analogous to the transverse instability 
of a density depletion in a conservative Gross-Pitaevskii 
equation [17] : the motion of grey solitons is inversely pro- 
portional to their depth, so modulation in the transverse 
direction forces different parts of the front to move with 
different velocities leading to vortex pair formation. This 
suggests that the several sources of such flows may con- 
tinuously generate a large number of vortices leading to 
a turbulent flow. Another possibility to create a turbu- 
lent flow is related to the formation of vortex lattice in 
a harmonic trapping potential due to an instability of a 
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FIG. 2: (color online) The evolution of the density of vor- 
tices as a function of time. For time t < 500 the pumping 
is nonuniform as discussed in the text. At time t = 500 
the nonuniformity of the pumping field is removed and the 
vortex density decays linearly initially, as the inset shows. 
As the density of vortices decreases, the system reaches the 
superfluid turbulence regime of well-separated vortices with 
a logarithmic decay [20]. These decay rates are in contrast 
with a power-decay rates of the order t~ 3 ^ in classical 2D 
viscous fluids and in the limit of the cGLE equation with zero 
dispersion [2"T] . 
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non-rotating solution [S]. By removing the circular sym- 
metry of either the trapping potential or pumping field it 
is possible to create a turbulent now of vortices instead 
of a regular vortex lattice (see the right panel of Fig. [lj . 

Numerical set-up. In order to engineer a turbulent for- 
mation and interaction of vortices we shall consider an in- 
homogeneous pump a(x) that can be obtained by passing 
the laser beam through a spatial phase (light) modulator. 
This will be even further simplified by assuming that only 
two laser intensities are allowed: the background with a 
superimposed set of almost periodical spots of a higher 
intensity, so that cv(x) = a everywhere except for x in- 
side the circles |x — aj| 2 < a where a(x) = ct\ + ao and 
Cj = c+Xi, &i = (ai±iT+Si, O2±iT+0j), T is the period 
of spots, c is the square of the spot radius, i = 0, 1, 2... 
and Xii $i an d <fii are random displacements of the order 
of the healing length. Both r\ and a take different val- 
ues for different pumping intensities. In practice, setting 
different values for these quantities does not change the 
qualitative behaviour of the system. In what follows both 
r\ and a with be set to be constants across the fields [19] . 

Through the time evolution one can observe the for- 
mation of vortices until their number begins to fluctuate 
about a constant value; see Fig. [2j The larger differ- 
ence between the two pumping intensities, a±, leads to 
the faster outflows and a larger number of vortices gener- 
ated. The relaxation has a negative effect on the number 
of vortices (compare the vortex densities for a\ = 3 and 
7} = or r) = 0.01 on Fig. At time t s = 500 (well after 



the quasi-equilibrium is reached) we remove the nonuni- 
formity of the pump by setting ct\ = 0. After that the 
vortices start annihilating each other leading to the decay 
of the turbulence. This stage can be compared and con- 
trasted with the wave turbulence of the Gross-Pitaevskii 
equation where the dissipation is at a given (high) mo- 
menta and so has a different physical meaning |20| . 

By tuning the nonuniformity of the pumping field it 
is possible to reach different turbulent regimes. If the 
difference between intensities, a,\, is below a threshold 
or the distance between the spots of higher intensity is 
large, no vortices will be created. In a case of a moderate 
ai and only few spots a set of several well-formed well- 
separated vortex pairs is created and the system is in a 
superfluid turbulence state (see the left panel of Fig. [I] 
and the left inset of Fig. [3]) . By increasing the difference 
between intensities a\ it is possible to create the state of 
strong turbulence (where vortex cores start to overlap; 
see the top inset of Fig. [3]). It is, therefore, tempting 
to see if the system can be driven even further to enter 
the regime of weak turbulence in which all coherence is 
lost and all Fourier amplitudes have random phases. To 
verify this we calculated the second moment of the cor- 
relation function g 2 = (|V ; | 4 )/(IV'| 2 ) 2 - By Wick's theorem 
the state of the weak turbulence corresponds to 172 = 2. 
As shown on Fig. [3] by raising a\ it is possible for the 
system to reach the weak turbulence state. Note that 
the relaxation r\ increases 172- This occurs because the 
relaxation increases the rate at which vortex pairs anni- 
hilate by bringing the vortex cores closer to each other; 
this effect can be seen on Fig. [2] showing the number of 
vortices in quasi-equilibrium. The energy released from 
vortex annihilation becomes converted into acoustic en- 
ergy therefore increasing 172 • 

In order to describe the turbulence in the Eq.Q 
we shall assume that there exists an inertial range in 
the momentum space and that the role of pumping 
and dissipation is insignificant there. The evolu- 
tion equation for the wave spectrum defined by 
( a k 1 0'k) = n^S (ki — k.2), with a k being the Fourier 
transform of ip and are discrete wave vectors, can 
be obtained by using a random phase approximation 
and expanding in small nonlinearity [22] . The equation 
takes the form d t n kl (t) = Jd 2 k 2 d 2 k 3 d 2 k4W kl .jc 2 ;k 3 M x 
(n k3 n k4 n kl + "k 3 «k 4 ™k 2 - n kl n k . 2 n ks - n kl n ka n k4 ) , 
where W^m-MM = j0yzH^-i + k 2 - k 3 - k 4 ) 
8(k\ + k\ — k\ — fc|). Two solutions of this evolution 
equation correspond to a thermodynamic equipartition 
of the total kinetic energy E = J k 2 n k dk, so that 
n<k ~ k~ 2 and to an equipartition of the total number 
of particles N — Jn^rfk, so that n k ~ const. These 
correspond to the two limits of the Rayleigh- Jeans 
distribution T j(k 2 + fi), where T is the temperature. 

We verified the existence of the inertial range in 
our simulations. Although the system is in a quasi- 
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FIG. 3: (color online) The second moment of the correla- FIG. 4: (color online) The wave spectrum log(nfc) vs log(fc) 

tion function Q2 as a function of ot\. Starting from an initial for the state of strong turbulence established for parameters 

constant density profile the nonuniform pumping is applied, n — 0, a — 0.3 at t = 450 (nonuniform pumping field with 

so that <?2 rapidly grows reaching a quasi-stationary state af- at\ — 3, ao = 1/2, top blue curve) and at t = 650 (uniform 

ter t ~ 20. After that <?2 fluctuates about a constant value. pumping filed ai = 0, ag — 1/2, bottom red curve). Lines 

The blue dots (n = 0) and red squares (77 = 0.01) show the corresponding to nt ~ const and rik ~ fc~ 2 are included, 

average of the g-2 during the time interval [20, t 3 ]- The time Inset shows the wave spectrum for the state of superfluid tur- 

snapshots of the normalized density l^ 2 of the fields if> are bulence ot\ — 1 with nu ~ k~ 2 spectrum of the inertial range 

shown for superfluid turbulence state (left inset), strong tur- and for the transitional state qi = 2, both for t <t 3 . 
bulence (bottom inset) and weak turbulence state (top inset) 




equilibrium rather than in the true thermodynamical 
equilibrium we observed both spectra. For the nonuni- 
form pumping the wave spectrum shows the particle 
equipartitions for strong turbulence (see the top (blue) 
curve of Fig. [2]), whereas the superfluid turbulence spec- 
tra corresponds to energy equipartition (the Kolmogorov- 
Zakharov energy cascade rik ~ k ~ 2 ); see the grey (green) 
curve on the inset of Fig. [I] During the turbulence de- 
cay stage the wave spectrum corresponds to rik ~ k~ 2 : 
see the bottom (red) curve of Fig. |2j This suggests that 
at these intermediate scales of the inertia range the tur- 
bulence is dominated by four-wave interactions and the 
wave field is weakly nonlinear and dominated by acoustic 
modes. 

In summary, we proposed a way to generate vari- 
ous regimes of turbulence in nonequilibrium condensates, 
such as exciton-polariton condensates. By designing a 
nonuniform pumping field that leads to sufficiently strong 
interacting fluxes it is possible to create the superfluid 
turbulence with well separated quantised vortices, the 
strong turbulence with overlapping and de-structured 
vortices or the weak turbulence state with a complete loss 
of coherence. The nonequilibrium condensates, therefore, 
are new and exciting systems with a nontrivial evolution 
of complex matter field with turbulence that may span 
regimes fundamentally different from the classical fluid 
turbulence. 

The author acknowledges useful discussions with 
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